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I. INTRODUCTION 



In this paper we conduct an investigation of the behavior of a periodically driven nonlinear 
mechanical system. The system is of interest both from the point of view of the relatively 
rich structure it possesses and also because it is fairly easy to construct a mechanical model 
of it. 

We begin with some standard techniques used to analyze such systems. Bifurcation dia- 
grams are plotted which locate regions of quasiperiodic, periodic and chaotic behavior within 
the parameter space of the system. Dimensions of the attractor are calculated for different 
values of control parameters and compared with dimensions estimated from Lyapunov ex- 
ponents. Numerical results obtained are in agreement with the Kaplan- Yorke conjecture. 
Multiple attractors are found in different regions in parameter space. The most interesting is 
the coexistence of a chaotic attractor and a limit cycle, which implies that initial conditions 
may determine whether the system will exhibit periodic or chaotic behavior. Boundaries of 
the basins of attraction appear to be smooth for high dissipation in the system, and they 
show fractal structure as the dissipation decreases. Creation and destruction of multiple 
attractors in this model is caused by tangent bifurcations and crisis phenomena. This illus- 
trates how even a simple nonlinear mechanical system may exhibit fairly complex behavior 
and posses a variety of chaotic features. 

We also conduct a symbolic analysis of the system. Symbolic dynamics provides almost 
the only rigorous approach to study the motion of dynamic systems. Such an analysis of 
one-dimensional maps on the interval is well understood Jl], 0. For the simplest case of the 
unimodal map, a binary generating partition may be introduced by splitting the interval 
at the critical point. By assigning the letter i? or L to a point of an orbit, depending 
whether it falls to the right or left side of the critical point, the orbit can be encoded with 
a symbolic sequence. The kneading sequence, which is the forward sequence of the critical 
value, determines all the admissible sequences of the map. By means of the kneading theory, 
all sequences are well ordered. The kneading sequence is the greatest. There is no allowed 
sequence which has a shifted sub-sequence greater than the kneading sequence. 

The symbolic dynamics of the unimodal map can be extended to that of two-dimensional 
(2D) maps. Two much-studied 2D models are the Henon map ||^ and its piecewise linear 
version, the Lozi map [|, ^. First of all, one needs to construct a 'good' binary parti- 
tion for them. In Ref. |^], by considering all 'primary' homoclinic tangencies of the Henon 
map, a method for determining a partition line was proposed. Once the binary parti- 
tion is determined, any orbit may be associated with a doubly infinite symbolic sequence 
S = ■ ■ ■ Sjn - ■ ■ S2Sj»SoSi ■ ■ ■ s„ ■ ■ ■ where Sq indicates the code of the initial point. The forward 
sequence •SqSi ■ ■ ■ s„ ■ ■ ■ and backward sequence ■ ■ ■ Sm - ■ ■ SgS-* correspond to the forward 
and backward orbits of the initial point, respectively. In Ref.0 the ordering rules for for- 
ward and backward sequences were discussed, and by introducing a metric representation 
of sequences the symbolic plane was constructed. It was pointed out that every primary 
homoclinic tangency cuts out a rectangle of forbidden sequences in the symbolic plane. The 
rectangles so deleted build up a pruning front which is monotonic across half the symbolic 
plane. By generalizing stable and unstable manifolds of an unstable fixed or periodic point 
to forward and backward foliations of any point 0, homoclinic tangencies are generalized 
to tangencies between the two classes of foliations. In terms of the generalization, we can 
make symbolic analysis of the orbits not only in but also out of attractor, including transient 
orbits. The extension of symbolic dynamics of maps from ID to 2D is made by decompos- 
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ing a 2D map into two ID maps based on forward and backward foliations. The coupling 
between the two ID maps is described by the pruning front or the symbolic representation 
of the partition line. An attractor is associated with backward foliations. 

It has been convincingly demonstrated that a properly constructed two-dimensional sym- 
bohc dynamics, being a coarse-grained description, provides a powerful tool to capture 
global, topological aspects of low-dimensional dissipative systems of ordinary differential 
equations (ODEs) g |TU|, [TT], |g, 0, 0, |T|, |g, [T^, |Tg. Therefore, to have a further 
global understanding of the bifurcation and chaos "spectrum" in the system we have studied 
i.e. the systematics of stable and unstable periodic orbits at varying and fixed param- 



eters, the types of chaotic attractors, etc., we shall carry out a full analysis of symbolic 
dynamics for this mechanical model. Based on the primary tangencies between forward and 
backward foliations in the Poincare section the appropriate partition of the phase space is 
made and three letters are used to describe the dynamics. From the ordering rules and met- 
ric representation of forward and backward sequences symbolic planes are constructed and 
the admissibility conditions for allowed sequences derived. Unstable allowed periodic orbits 
embedded in a chaotic attractor are predicted through symbolic analysis and verified nu- 
merically as well. In the parameter space period windows are located and the corresponding 
period words are determined according to the partitioning of phase portraits. 

The paper is organized as follows. In section Section II a brief description of the me- 
chanical system and a derivation of the equations of motion is given. In Section |T| chaos is 
detected, and regions of periodic, quasiperiodic and chaotic behavior are located within the 
parameter space. Lyapunov exponents and associated dimensions are calculated for certain 
attractors, and the Kaplan- Yorke conjecture is established. In section |IV| crisis phenom- 
ena and coexisting attractors are studied. We next turn to the symbolic dynamics of the 
module. In Section |V|, the global dynamical behavior of the system while parameters vary 
are shown in phase space, and the partition lines are determined from tangencies between 
forward and backward foliations in the Poincare section. In Section by introducing a 
metric representation of forward and backward sequences according to their ordering rule 
the symbolic plane is constructed, and the admissibility condition of sequences is discussed. 
We then analyze the admissibility of periodic sequences based on a finite number of points 
on the partition lines, and give the numerical results of allowed periodic orbits in Section 



VI I| . We also locate the period windows for some range of parameters, assign words to the 
stable period orbits based on the partitioning of phase space, and discuss their ordering 
properties as parameters change in Section |VIII| . In Section |^ we discuss some aspects of 
this model as used in a mechanical propulsion device. Finally, in Section ^ we present some 
conclusions. 
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II. THE EQUATION OF MOTION 




FIG. 1: Two gears and the rod 



The system we study consists of two gears and a rod. The first one, the solar gear, has 
radius i?', and is fixed so neither its center of mass can move nor can it rotate around its 
axis. The second, the planetary gear, is of radius R and is attached to the first gear so when 
it rotates around its axis at the same time it goes around the fixed gear. Wc assume that 
the movable gear is powered by some device which keeps it moving with constant angular 
velocity ujd- When the planetary gear rotates for the angle (p around its axis it moves around 
solar gear for an angle (f)' (Fig. II). The relation between two angles is 

= R(j) (1) 

A uniformly dense rod of length L is at one end joined to the fringe of movable gear. The 
rod can rotate around the joint, assumed without friction. Both gears and the rod lie in a 
horizontal plane so the potential energy of the system is constant during the motion. We 
choose a coordinate system with origin at the center of the fixed gear. The system has two 
degrees of freedom, and a convenient choice of coordinates are 0', the angle that the second 
gear makes with x-axis, and i9, the angle that the rod makes with the x-axis. Since the 
potential energy for the system is constant as well as the angular velocity of the planetary 
gear, the Lagrangian for the system is determined by the kinetic energy of the rod only. A 
general expression for the kinetic energy for the rigid body that moves in two dimensions is 

T^^Jdmix + y) (2) 

We may express the position coordinates for the mass element dm over coordinates 4> and 
■J? as 

x^ {R' + R) cos (p' - R cos(</)' + 0) + / cos , . 
y^{R' + R) sin (p' - R sin(0' + 0) + Z sin i? 
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where / is position of the mass element dm at the rod. Deriving x and y over time and 
ehminating 0' rather than 0, using ([^), we get: 

X = —(1 + r)R(j) sin r(f) + R{1 + r)0sin [(1 + r) 0] — Wsiwd , ■. 

y = {1 + r)_R0cosr0 — i?0 (1 + r) cos [(1 + r) 0] + Wcos-t} 

where r = R/R'. Our initial assumption was that angular velocity of the planetary gear 
was constant = uod. For simplicity, we will assume that initial condition is chosen so that 
= ujdt. Substituting derivatives @) into (H) and integrating over the whole rod, assuming 
its linear mass is constant, we get the kinetic energy, and therefore the Lagrangian, to be: 

1-1 

L = T= + -(1 + r)mLRud^ {cos{^d - rcj)) - cos - (I + r) 0]} (5) 

Here / is the moment of inertia of the rod, and we dropped terms that can be written as 
a total time derivative of some function and therefore do not contribute to the Lagrangian. 
Substituting this Lagrangian into the general expression for the Euler-Lagrange equations 
of motion of dissipative systems 

d dL dL 

and choosing and for coordinates we obtain the equation of motion to be 
1 + r 

M H —mLRjjj {r sin(t9 - r0) - (1 + r) sin [?9 - (1 + r)0]} = (7) 

where is some friction force that acts upon the rod, and / is its moment of inertia. If 
we assume that the friction, as is commonly done for a pendulum, has the simple form 
= —rjd (r/ > 0), we may write the equations of motion as 

d + + uo^ |sin(?9 - r0) - sin[?9 - (1 + r)0]| = 0, (8) 

where 

^o=r{l + r)^ul (9) 

If we introduce now a new variable 9 = d — rcf) we can bring the equations of motion to the 
arguably simpler form: 

e + ^e + sin e = -'^ujd + ^^c^n sin(^ - 0). (10) 
/ I r 

This equation is nonlinear with a coupled term, and must be solved numerically. Therefore, 
for convenience we will write it in a dimensionless form by introducing the dimensionless 
time variable r = ujot. Equation (|iy) then becomes 

^ + -^ + sin^ = -— + sin(^-0), (11) 

Q Q r 



where 



Ud / 2L 



LUo V3r(l + r)i? 
Q = ^f^^MI, (12) 



<x> 

■a 
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FIG. 2: (a) Bifurcation graph showing broad chaotic region over segment of Q. a = 1.2, r = 0.9. 
(b) Corresponding change of the largest Lyapunov exponent. 



and 

= LJdt = ar. (13) 

Since we assume throughout this paper that ud is constant, the state of this system can be 
uniquely described in 3-dimensional phase space {6, 6, 0). Furthermore, from the fact that 
depends hnearly on time we may infer that trajectories in phase space will be "smooth" in 
the direction, and that all chaotic features, if any, can be observed in the Poincare section 
= constant. In this way the problem of analyzing this three dimensional dynamical flow 
is reduced to the analysis of a two dimensional map. 

Equation ( pT]) has a similar form to that of the equation of motion for the driven pendu- 
lum, as well as to the Mathieu equation; however, there are some essential differences. The 
virtual drive frequency a in the normalized equation (|ll]) does not depend upon the drive 
frequency - the angular velocity of the moving gear - as it is constant which depends 
only on the dimensions of the rod and gears. On the other hand, the "quality" factor Q 
is proportional to cj^, so a change of drive frequency will eventually appear as a change of 
dissipation in the system. 



III. ROUTES TO CHAOS 

For our system, we observe three characteristic regions in the parameter space - quasiperi- 
odic, periodic and chaotic. In order to locate these regions we investigate the parameter 
space using bifurcation diagrams. Such diagrams show the evolution of the attractor for a 
dissipative system with a change of the system's parameters. We choose to plot a projection 
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FIG. 3: Power spectra for different values of quality factor Q: (a) Q = 0.2 motion is quasiperiodic, 
(b) Q = 0.9, motion is periodic with period one, (c) Q = 1.8 and (d) Q = 6.5, motion is chaotic. 
G(l(j) is the Fourier transform of 0{t). The frequency on the x-axis is expressed in terms of the 
virtual drive frequency a. 



of the Poincare section (/) = on the 9 axis against the system's parameter Q, while keeping 
other two parameters fixed. We repeat this procedure for various values of a and r. 

With the increase of Q from zero on, the motion of the system is quasiperiodic, which is 
observed as a smeared region in the bifurcation graph. For certain value of Q the motion 
suddenly becomes periodic, and with a further increase of the quality factor Q the system 
reaches chaos through a sequence of period doublings. For the parameters a,r < 0.5, the 
quasiperiodic region is relatively short, and chaos is found only in a narrow region of Q 
values. The richest chaotic structure is found for the choice of parameters a, r ~ 1. Finally, 
for a, r > 1 quasiperiodic motion becomes dominant over the range of Q, and periodic and 
possible chaotic behavior are found only in narrow isolated regions. 

Solid proof of chaos is the existence of at least one positive Lyapunov exponent. The 
spectrum of Lyapunov exponents gives us, besides qualitative information about a system's 
behavior, also a quantitative measure of the system's stability [^, |2^, Evolution 



of the Lyapunov exponent spectrum with a change of parameters gives us additional infor- 
mation about the system's behavior to that from bifurcation graphs. Combined, these two 
methods are reliable tools for examining a system's behavior and routes to chaos. Figure |^ 
gives a comparison of the bifurcation diagram and Lyapunov spectrum evolution for a range 
of Q. Since the system is represented in 3 dimensional space it is sufficient to determine 
only the largest Lyapunov exponent. The other two can be inferred from Haken's theorem 
p5| , and the fact that the dissipation of the system is constant throughout phase space, so 
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FIG. 4: (a) Poincare section of the quasiper iodic attractor is a smooth hne. a = 1.2, r = 0.9, 
Q = 0.5. (b) Poincare section of the chaotic attractor reveals a kneading sequence. Q = 1.8. (c) 
Higher dimensional chaotic attractor. Q = 6.5. 



the sum of Lyapunov exponents must equal to its numerical value J2i = ^/Q- 

A detailed description of the parameter space for this system is given in Ref. [^. Here 



we illustrate its main features with a bifurcation diagram plotted for parameters a = 1.2 
and 0.9 (Fig. 0a). The smeared region in the bifurcation diagram that is observed for the 
values < Q < 0.6378 . . . corresponds to a quasiperiodic behavior of the system. Although 
it is generally possible to notice differences between chaotic and quasiperiodic regions in a 
bifurcation graph, it is necessary to employ other methods in order to have solid evidence 
of a system's behavior. 

Analysis of the power spectrum of the numerical solution is one of the methods that is 
particularly suitable for examining quasiperiodic motion. A power spectrum for Q = 0.5 
(Fig. ^) shows that the motion of the system is a superposition of a finite number of 
harmonic modes, and therefore is not chaotic. The motion is not periodic though, since fre- 
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Q 


Dl 


Do 


Dl 


D2 


0.5 


2.000 ± 0.001 


2.002 ± 0.003 


2.007 ±0.006 


1.999 ± 0.002 


1.8 


2.14 ±0.01 


2.11 ±0.01 


2.044 ± 0.007 


2.012 ±0.003 


6.5 


2.41 ±0.01 


2.39 ±0.01 


2.389 ± 0.004 


2.299 ± 0.001 



TABLE I: Various dimensions of the attractors. 



quencies of harmonics are incommensurate, and the trajectory in phase space never retracts 
itself. The attractor for the system is a smooth 2-dimensional surface in phase space. The 
Poincare section for Q = 0.5 (Fig. fTTja) indeed reveals smooth attractor. 

For values of 0.7442 . . . < Q < 0.9300 ... it is clear from the bifurcation graph that the 
system's behavior is periodic with period that is equal to the virtual drive period 271 /a. This 
is also seen in the power spectrum for Q = 0.9 in Figure ^d. With an increase of Q over 
1.2040. . . the system becomes chaotic through a period doubling sequence. The smeared 
region in the bifurcation diagram indicates that the attractor for the system has a complex 
structure in that region of the parameter space. The power spectrum also shows complexity 
of the system's behavior and indicates chaos in that region of parameter space. The power 
spectrum does not consist of a finite number of harmonics and their integer multiples any 
more, but is rather a broad band over the frequency axis (Figs. ^ and |^). In order 
to obtain a solid evidence of chaos one has to show that at least one Lyapunov exponent 



for the system is positive in that region |27|. Indeed, Figure 0b shows that the largest 



Lyapunov exponent is positive in this region of parameter space, apart from windows of 
periodic behavior. 

In particular, at Q = 1.8 and Q = 6.5 the largest Lyapunov exponents are found to be 
Ai = 0.093 ± 0.005 and Ai = 0.106 ± 0.005, respectively. The system's attractors for these 
values of parameters (Figs. |ITJ) and |IT|c) show characteristic stretching and folding pattern, 
indicating their fractal structure. 

In the quasiperiodic region the attractor has dimension 2, as the motion is performed 
over a smooth surface, while the limit cycles have dimension 1. Since motion takes place in 
3-dimensional phase space it is reasonable to expect that strange attractors have dimension 
between 2 and 3. We can use the fact that the attractor is always smooth in the direction, 
so in order to estimate its dimension we can estimate the dimension of its Poincare section 
and obtain the attractor's dimension simply by adding one to it. 

In order to estimate the dimension of an attractor lying in n-dimensional phase space, the 
attractor, first, has to be covered by a grid of ?7,-dimensional hypercubes of size e, and then 
the probability of finding a point of the attractor in each hypercube Pj has to be determined. 
The index i here refers to a particular cube. In our case we consider a 2-dimensional Poincare 
section, so the hypercubes are actually squares. A general expression for the dimensions of 
the g*^-order is then given by: 

1 In P '' 

^. = lim 7^^, (14) 

e^o q — I In e 



where the summation is over all hypercubes where Pj > ||28[ . The parameter q ranges from 
—oo<q<oo, and for qi > q2 we have Dg^ < Dg^. The most commonly used dimensions 
are Dq, Di, and D2, or the capacity, information, and correlation dimension, respectively. 
We will limit ourselves to estimating these three dimensions only. 
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On the other hand, the Lyapunov dimension of the attractor is defined as: 



D,=3 + + + (15) 

where j is a largest integer for which Ai + A2 + • • • + Aj > 0. Kaplan and Yorke conjectured 
that the Lyapunov dimension represents an upper limit to the information dimension, i.e. 
Dl> Di It is to be expected though that numerical values of these four dimensions 

fall pretty close to each other 



We test this conjecture on our calculations for the quasiperiodic attractor (Fig. |TT|), and 
as expected we obtain the Lyapunov, capacity, information, and correlation dimensions to 
be equal 2, within the limits of uncertainty of our numerical methods. Furthermore, we 



can estimate these dimensions for the attractors in Figures |T|b and |T|c. Both chaotic 
attractors have noninteger dimensions (Table H), confirming their fractal structure. In all 
three cases the inequality Dl > Di holds, suggesting that this system conforms Kaplan- 
Yorke conjecture. 



IV. COEXISTING ATTRACTORS AND CRISIS 



e 
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FIG. 5: Bifurcation diagram indicating coexisting attractors. Dashed line shows unstable limit 
cycle. 



An interesting feature of this system is that it often has periodic and chaotic attractors 
coexisting at the same point in parameter space. This means that the stability of this 
system may depend on the initial conditions - whether we choose it in the basin of periodic 
or chaotic attractor. Here we shall present two examples which illustrate the importance of 
understanding system's dynamics in order to properly assess its stability. 

The bifurcation diagram |^ shows that a tangent bifurcation ||32| occurs Q = 1.0105 . . ., 
which creates a pair of stable and unstable limit cycles in the phase space, while a chaotic 
attractor still exists. Accordingly, the phase space splits into a basin of the stable limit 
cycle and a basin of the chaotic attractor. The boundary between the two basins appears 
to be a smooth curve (Fig. With an increase of the control parameter Q the basin 
of attraction for the limit cycle expands until the unstable limit cycle, which lies on its 
boundary, collides with the chaotic attractor. The chaotic attractor experiences a boundary 
crisis [0 and disappears along with its own basin. The basin of the limit cycle then suddenly 
expands, occupying the whole phase space. With a further increase of the control parameter 
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-3 3 

9 

FIG. 6: Poincare section of the chaotic attractor at the verge of crisis, shown within its basin of 
attraction in. White dots represent a stable, and black dots an unstable, limit cycle. Q = 1.8. 



the limit cycle evolves into a chaotic attractor through a sequence of period doublings. 
The newly created chaotic attractor expands in size as the parameter Q increases, and 
eventually collides with the unstable limit cycle. The system comes to a crisis again, but 
since the unstable orbit is not at the basin boundary any more the crisis is internal, and the 
attractor suddenly expands its size [Q. This expanded attractor contains features of the 
chaotic attractor destroyed in the boundary crisis, indicating that information about the 
system's behavior at lower values of Q has "survived" the crisis. 



-a 
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FIG. 7: Bifurcation diagram shows creation of coexisting attractors. 
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FIG. 8: Three coexisting attractors. Period-three hmit cycle (black points) and period-one limit 
cycle (white point) are enlarged for clarity. Q = 3.515. 



As we increase the control parameter Q further, we observe another occurrence of coex- 
isting attractors (Fig. 0). Tangent bifurcations aX Q = 3.2398 . . . and Q = 3.4177. . . create 
stable-unstable pairs of period-one and period-three orbits, respectively, so we have three 
coexisting attractors. Basin boundaries are now locally disconnected fractal curves |33], as 



shown in Figure |^. Basins of attraction are highly interwoven in this region, and for prac- 
tical purposes it is difficult to determine initial conditions that would lead to a particular 
attractor. The newly created period-3 orbit quickly evolves to a chaotic attractor with an 
increase of the control parameter (Fig. |^), and disappears in a boundary crisis. Basins of 
attractions in Figures 10 and p| are plotted using software package Dynamics 2 BH] . 



V. PARTITIONING OF THE POINCARE SECTION 

We now turn to a symbolic analysis of this model. For this, we shall choose convenient 
phase space x = 6/27i, y = x, and write the equations of motion in their nonautonomous 
form as 

x = y, /^gN 

y = -y/Q +[-ar/Q -sm{2TTx) + {l + r)/rsm{27TX - aT)]/{2'K). ^ ' 

To show its richer properties in the phase and parameter spaces we fix r = 1.088 and vary the 
values of Q and a in such a way that a = 0.8 + 0.3Q. All Poincare sections are constructed 
in (j) = 277r/64 plane of phase space. 

To observe the dynamics changes for varying Q and a we juxtapose the Poincare sections 
(x, y) of attractors, the Xn — > Xn+i ffist return map and the fi — u portraits where n = 
1 - (0.4 + y) cos(27ra;) and z/ = 1 - (0.4 + y) sin(27rx) at Q = 0.57, 0.76, 1.2577, 1.68, and 2.0 
from top to bottom in Fig. ^. Obviously, when Q increases, 1) the portraits undergo the 
changes: ID ^ 2D — > ID; 2) the return map Xn+i — Xn shows: a subcritical circle map 
without any decreasing branch — > the appearance of decreasing branches — > a subcritical 
circle map; 3) the dynamical behavior exhibits: quasiperiodicity chaos quasiperiodicity. 

The quasiperiodic attractor and forward foliations aX Q = 0.57 are shown in Fig. [1^, 
from which one can see clearly that there is no any tangency between the attractor (part of 
backward foliations) and the forward foliations. As a matter of fact, when Q increases to a 
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FIG. 9: Attractors shown in the x — y, Xn — Xn+i, and ^ — v planes at 
0.57, 0.76, 1.2577, 1.68, and 2.0 from top to bot^c^n for r = 1.088, and a = 0.8 + 0.3Q. 




FIG. 10: Attractor (dots) and forward foliations (dash curves) for Q = 0.57. 



value the tangencies begin to appear, and correspondingly one can obtain a critical circle or 
annular map Xn+i — Xn- If Q keeps increasing, then the map would be supercritical. Finally 
there are no tangency and the map becomes subcritical again (see Fig. P). This shows the 
close connection between the geometric properties (the tangencies) of an attractor and the 
dynamical behavior (quasiperiodicity and chaos) of a system. 

In Fig. [Tl| we show the attractor and two primary partition lines {•B and •C) on the 
background of the forward foliations (dash curves) for Q = 0.76. The line marked with •A 
is the pre-image of •B. The areas in between these lines are labeled by •R, •L, and •N. The 



corresponding first return map is plotted in Fig. |T2|. The three monotone segments in Fig. ^ 
are assigned the letters L, R, and N, in accordance with the two-dimensional partitions in 



Fig. 11 



A more interesting case is encountered a.t Q = 1.2577. The attractor, three forward 



foliations passing tangencies and the primary partition lines are shown in Fig. O, which 



manifestly exhibits two-dimensional feature. We also separately plot in Figs. and [T^ the 
Xn+i — Xn first return map constructed from Fig. |T3| by using the x coordinates and its 
swapped form x'^^^ — with 



X 



x + 1 — xb for X < xb, or x' 



X 



Xb for X > Xb, 



(17) 



where is the x coordinate of the attractor tangency on the partition line •B in Fig. [13. 
Actually the partition lines and •A are not respectively parallel to the x^+i-axis and 



x„-axis in Figs. |IJ and [15. 
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FIG. 11: The Poincarc map (dots) and forward foliations (dash curves) at Q = 0.76. The primary 
partition Hnes •B, and the pre-image •A of divide the attractor into three parts labeled by 
the letters L, R and N. 



VI. ORDERING RULE AND ADMISSIBILITY CONDITION 

After the partition lines are determined, each point, or its orbit, may be encoded 
with a doubly infinite symbolic sequence consisting of the letters R, L and N, say, 
S — ■ ■ ■ Sjn - ■ ■ S2S1 •S1S2 • • • Sn • • where s„ is the code for the n-th point of the forward orbit, 
and the code for the m-th point of the backward orbit. The "present" position is indicated 
by a solid dot, which divides the doubly infinite sequence into two semi- infinite sequences, 
i.e., the backward sequence ■ ■ ■ Sm - ■ ■ Sf-^i* forward sequence •S1S2 ■ ■ ■ s„ ■ ■ ■. 

A metric representation for symbolic sequences can be introduced by assigning numbers 
in [0, 1] to forward and backward sequences. We first assign an integer = —1 or 1 to the 
symbol Si when it is the letter N or otherwise. Then we assign to the forward sequence 
•S1S2 • • • ■ ■ ■ the number 

CO 

a = J2f^i^~\ (18) 

i=l 

where ^ 

l^i=\l ior Si^{R if n = 1' (19) 



or, 

l_ii = {l ior Si={R if n = -1- (20) 









if lie, 
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if n^. 
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FIG. 12: The Xfi^i — Xfi first return map at Q — 0.76. 
Similarly, the (3 assigned to the backward sequence • • • Sm • • • S2Sj» is defined by 



where 



or 





/5 = E^3-^ 






i=l 






(R 






for Sj = < AT 


and n = 1' 




[l 






(R 


i-l 




for < N 


and n = ~1 


lo 


[l 





According to the definition we have 

a{»NL°°) = p{L°°m) = 1, a{»L°°) = (5{R°°») = 0, 

a{»NNL°°) = a{»RNL°°) = 2/3, P{R°°L») = /3(i?~A^«) = 2/3, 

a{»RL^) = a{»LNL°°) = 1/3, P{L°°N») = P{L^Rm) = 1/3. 



(21) 



(22) 



(23) 



(24) 



In this representation a bi-infinite symbolic sequence with the present dot specified corre- 
sponds to a point in the unit square of the a — (3 plane, the so-called symbolic plane. In the 
plane, forward and backward foliations become vertical and horizontal lines, respectively. 
We may define the ordering rules of forward (or backward) sequences according to their a 
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FIG. 13: The chaotic attractor (dots), forward fohations (dashed hnes), and partition hnes "A, 
•B, and #0 at Q = 1.2577. 

(or f3) values. From Eqs. ([T8|-8) we then have 

•EL--- <»ER--- <»EN---, •0L---> •0R---> •ON (25) 

and 

---RE^<---NE^<---LE^, ---R0^> ---N0^> - --LO^, (26) 

where the finite strings E and O consist of letters L, R and and contain an even and 
odd number of letter N, respectively. This ordering rule is similar to that for sequences of 
the dissipative standard map at some values of the parameters |]lOl . 



When foliations are well ordered, the geometry of a tangency places a restriction on 
allowed symbolic sequences. A point on the partition line C« (image of •C) may symbolically 
be represented as QC • P. The rectangle enclosed by the lines QN^, QR^, •P, and •NL°° 
forms a forbidden zone (FZ) in the symbolic plane. Therefore, a symbolic sequence IJ 
with /• between QN^ and QR^, and at the same time •J > •P must be forbidden by the 
tangency QC • P. In the symbolic plane the sequence IJ corresponds to a point inside the 
forbidden zone of QC • P. Similarly, UB •V stands for a tangency on the partition line 
(image of •B). The lines UN^, •¥ and •L'^ enclose a rectangle FZ in the symbolic 

plane. Any sequence KT with between UL^ and UN^ while •T < •V is forbidden by 
the UB • V . Each tangency point on a partition line rules out a rectangle in the symbolic 
plane. The union of the forbidden rectangles, determined from one and the same partition 
line, forms the fundamental forbidden zone (FFZ), a boundary of which is the so-called 
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^n+1 




FIG. 14: The Xn+i — Xn first return map constructed from Fig. |13|. 



pruning front. Consider a finite set of tangencies {QiC • Pi} (or {UjB • Vj}). If the shift 
of a sequence ■ ■ ■ Sk-2Sk-i • SkSk+i--- satisfies the condition that the backward sequence 
• • ■Sk-2Sk-i* is not between QiN* and QiR* (or UjN* and UjL») , and at the same time 
•Pi > •SfcSfc+i ■ ■ ■ (or rnVj < •SkSk+i ■ ■ ■) for some i (j), then this shift is not forbidden by 
any tangencies of C« or B», owing to the property of well-ordering of foliations. Thus, we 
may say that the shift is allowed according to that tangency. A necessary and sufficient 
condition for a sequence to be allowed is that all of its shifts are allowed according to the 
two sets of tangencies. To check the admissibility condition, we consider again the two cases 
Q = 0.76 and Q = 1.2577, and draw 10000 points representing real sequences generated 



from the Poincare map together with the FFZ in the symbolic plane, as seen in Figs. |T6| and 
T7|. One can see that the FFZ indeed contains no point of allowed sequences. A blow-up 



of the right-hand side pruning front in Fig. |1^ is displayed in Fig. |T8[ The structure means 
two-dimensional feature, related to the two tangent points in the upper part of attractor on 
the partition line (see Fig. |1^). We shall use the two tangencies to make 2D analysis of 
periodic sequences later on (see T3 and T4 in the next section). 



VII. UNSTABLE PERIODIC ORBIT SEQUENCES 



The attractor at Q = 0.76 does not show much two-dimensional nature, so the reduction 
to symbolic dynamics of one-dimensional circle map may capture much of the essentials. 
We start from this simple case. The attractor resembles that of a one-dimensional circle 
map except for a segment with two sheets, one of which without the part ( see Figs. O 
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FIG. 15: The x' 



swapped return map constructed from Fig. 14 



and |12D. From the two primary partition hnes •C and •B in Fig. we get the following 
sequences for attractor points: 

Ti : ■ ■ ■ A^A^ ■ ■ ■ NNAB • LRRLNLNLNLRRRLLRRR ■ • • , 
T2 : ■■■NN--- NNLC • RRLRLN LRRRLRLN RLN L 

In order to reduce the two-dimensional attractor to one-dimensional return map, we need 
to determine two kneading sequences Kb and Kc- They are the forward sequences of Ti 
and T2, respectively. 



Kb 

Kc 



LRRLNLNLNLRRRLLRRR 
RRLRLN LRRRLRLN RLN L 



(27) 



Compared with the original 2D map, the ID circle map given by these Kb and Kc puts less 
constraints on allowed orbits. Since the attractor has only one sheet crossing each primary 
partition line nearly no difference between the ID and 2D maps can be recognized if the 
sequences of short periodic orbits are concerned. 

The knowledge of the two kneading sequences ( P7| ) determines everything in the symbolic 
dynamics of the circle map For example, one may define a rotation number W, also 
called a winding number, for a symbolic sequence by counting the weight of letters R and 
N in the total number n of all letters: 



W = lim -(Number of R and N). 

n— >oo fi 



(28) 
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FIG. 16: The symbolic plane at Q = 0.76. 10000 points of real orbits generated from the Poincare 
map are also shown together with the FFZ in which no point falls. 



Chaotic regime is associated with the existence of a rotation interval, a closed interval in the 
parameter plane 
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Within a rotation interval there must be well-ordered orbits. We can 
construct some of these well-ordered sequences explicitly, knowing the kneading sequences 
Kb and Kc. 

In our case it can be verified that the ordered periodic orbits {RL)°° and {RRL)°° are 
admissible. These two sequences have rotation numbers 1/2 and 2/3, so the rotation interval 
of the circle map contains [1/2, 2/3], inside which there are rational rotation numbers 3/5, 
4/7 and 5/8 with denominators up to 8. Their corresponding ordered orbits are {R^ LRL)°° , 
[R^ L{RL)'^]°° and [{R^ L)"^ RL]°° . A very easy way to construct a longer well-ordered peri- 
odic sequence with a given rational rotation number from two shorter well-ordered periodic 

Take, for example W 



sequences can be found in Fig. 12.9 of Ref. 



3/5, 



3 
5 



2 
3 



2 + 1 

3 + 2' 



RRLRL = RRL + RL. 
We can further construct not-well-ordered sequences from well-ordered ones by the fol- 



lowing transformation. One notes that in Fig. |T2| the lower limit of •A is the greatest point 
on the subinterval L, while the upper limit of •A is the smallest R. When •A is crossed by 
a continuous change of initial points the corresponding symbolic sequences must change as 
follows: 

greatest LN ■ ■ ■ ^ smallest RL ■ ■ ■ . 
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FIG. 17: The symbolic plane at Q = 1.2577. Together with the FFZ, 10000 points representing 
real orbits are drawn. None of them falls inside the FFZ. 



Similarly, on crossing mC another change of symbols takes place: 

greatest R ^ smallest A^. 

Neither change has any effect on rotation numbers. As an example, starting with the ordered 
period 7 orbit [R'^{LRfL]°° we obtain 

RRLRLRL NRLRLRL NRLRLNL NRLNLNL NRLNRLL 

NRRLRLL RRRLRLL RRRLLNL ^ NRRLLNL NNRLLNL 

and 

NRLRLNL NRLRRLL RRLRRLL 

as candidates for the fundamental strings in not-well-ordered sequences of period 7. Among 
these sequences, {N N RLLN L)°° and {RRLRRLL)°° are forbidden by Kq and Kb, respec- 
tively. 

In this way we have determined all periodic sequences up to period 8, allowed by the 
two kneading sequences (|27|) . The result is summarized in Table ||. We have examined the 
admissibility of all these sequences by checking if their all shifts fall into the FFZ in the 
symbolic plane of Fig. |16[ They turn out to be totally allowed. In fact, by determining the 
symbolic sequence of every point in the attractor, we have numerically found all these orbits 
easily and listed the coordinate of the first letter in a sequence in Table 0. 
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0.65 



FFZ 





0.8794 



0.87965 



a 



FIG. 18: A blow-up of the symbolic plane Fig. 17 in the intervals a = [0.8794,0.87965] and 
P = [0,0.65]. 
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TABLE II: Allowed unstable periods up to 8 for Q = 0.76. The {x, y) is the coordinate of the first 
letter in a sequence. Only non-repeating strings of the sequences are given. P denotes the period 
and W the rotation number. 



P W Sequence 




X 




y 


2 


1/2 


RL 


0. 


,680655071 


-0, 


.184044078 


2 


1/2 


NL 


0. 


,269335571 


-0, 


.128679398 


3 


2/3 


RLR 


0. 


,666418156 


-0, 


.163903964 


3 


2/3 


RLN 


0. 


,679742958 


-0, 


.172914732 


4 


2/4 


NRLL 


0. 


,181303156 


-0, 


.181685939 


5 


3/5 


RRLRL 


0. 


,754061730 


-0, 


.232614919 


5 


3/5 


NRLRL 


0. 


,153005251 


-0, 


.198571214 


5 


3/5 


NRLNL 


0. 


,069305032 


-0, 


.242724621 


5 


3/5 


NRRLL 


0. 


,004459397 


-0, 


.266933087 


5 


3/5 


RRRLL 


0. 


,818461882 


-0, 


.262557048 


5 


3/5 


RRLNL 


0. 


,805821652 


-0, 


.257789205 


6 


3/6 


NRLRLL 


0. 


,162899971 


-0, 


.192734200 


6 


3/6 


NRLLNL 


0. 


,176235170 


-0, 


.184748719 


6 


4/6 


RRRLRL 


0. 


,876708867 


-0, 


.277086130 


6 


4/6 


NRRLRL 


0, 


,941828668 


-0, 


.278808659 


6 


4/6 


RRLNRL 


0. 


,781248148 


-0, 


.246935351 


7 


4/7 


RRLRLRL 


0. 


,753419555 


-0, 


.232248988 


7 


4/7 


NRLRLRL 


0. 


,159195421 


-0, 


.194929335 


7 


4/7 


RRLLNRL 


0. 


,744687447 


-0, 


.227153254 


7 


4/7 


NRLRLNL 


0. 


,151174731 


-0, 


.199641399 


7 


4/7 


NRLLNRL 


0, 


,185646405 


-0, 


.179051255 


7 


4/7 


NLNLNRL 


0. 


,259007234 


-0, 


.134664855 


7 


4/7 


RRRLLNL 


0. 


,850218105 


-0, 


.272009133 


7 


4/7 


NRRLLNL 


0. 


,981762965 


-0, 


.272665749 


7 


4/7 


RRRLRLL 


0. 


,858495739 


-0, 


.273871030 


7 


4/7 


NRRLRLL 


0. 


,964283072 


-0, 


.275994309 


8 


4/8 


NRLRLRLL 


0. 


,160408189 


-0, 


.194212064 


8 


4/8 


NRLRLLNL 


0. 


,162092260 


-0, 


.193213633 


8 


4/8 


NRLLNLNL 


0. 


,177466432 


-0, 


.184005818 


8 


5/8 


RRLRRLRL 


0. 


,759627678 


-0, 


.235734676 


8 


5/8 


NRLRRLRL 


0. 


,143963347 


-0, 


.203824391 


8 


5/8 


NRLNRLRL 


0. 


,093629167 


-0, 


.231107382 


8 


5/8 


NRRLRLRL 


0. 


,953347043 


-0, 


.277574465 


8 


5/8 


RRRLRLRL 


0. 


,869038801 


-0, 


.275879784 


8 


5/8 


NRLNLNRL 0. 


,070703636 


-0, 


.242089800 


8 


5/8 


NRLNRRLL 


0. 


,111829461 


-0, 


.221687365 


8 


5/8 


NRLRRRLL 


0. 


,131212194 


-0, 


.211077481 


8 


5/8 


RRLNLNRL 


0, 


,792783112 


-0, 


.252286399 


8 


5/8 


RRLNRLRL 




,782318121 


-0, 


.247450423 


8 


5/8 


NLNRRLRL 


0.351020175 


-0, 


.087984499 


8 


5/8 


NLRRRLRL 


0.353562509 


-0, 


.086982096 



For the more interesting case Q = 1.2577, based on Fig. ^ we list the following five 
tangencies along the B» and C« lines: 



Ti 


■■■L-- 


■ LLAB • 


LRNNRLLNNNNNRLRLLN ■ ■ 


T2 


■■■N- 


■ ■ NNRC 


• NRLLRNNNNRLRLRNNRN 


T3 


■■■L-- 


■ LLNC • 


NRLLNNNRNRNNRNRNNR 


T4 


■■■N- 


■ ■ NNNC 


• NRLLNNNRLRNNNNRLLN 


T5 


■■■N- 


■ ■ NRLC • NNNNRLRLNRNRLNNNRN 



From Ti and T2 whose forward sequence is the greatest among the tangencies along C«, we 

get 

Kb = LRNNRLLNNNNNRLRLLN 
Kc = NRLLRNNNNRLRLRNNRN ■■■ . 
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TABLE III: Allowed unstable periods up to 7 for Q = 1.2577 in the ID case; those with an asterisk 
are forbidden by 2D tangency T5. 



P W Sequences 

1 I/IR N 

2 2/2 RN 

2 1/2 RL NL 

3 3/3 NRR RNN 

3 2/3 RLR RLN 

4 4/4 RRRN RRNN RNNN 
4 2/4 NRLL 

4 3/4 RRRL NRRL NNRL* 

5 5/5 RRRRN RRRNN RRNNN RNNNN RNNRN RRNRN 
5 3/5 RRLRL NRLRL NRLNL NRRLL NNRLL* 

5 4/5 RRRRL NRRRL NNRRL NNNRL RNRRL RNRLN RNRLR 

6 6/6 RRRRRN RRRRNN RRRNNN RRNNNN RNNNNN RRRNRN 
6 6/6 RRNNRN RNNNRN NNRRNR 

6 3/6 RLRLLN RLLNLN 

6 4/6 RRRLRL RLNRRL RLNNRL* NLNNRL* RLRRLN RRLNLN 
6 5/6 RRRRRL N RRRRL N NRRRL NNNRRL NNNNRL* RNRRRL 

6 5/6 N RNRRL RRNRRL NNRNRL* RNRNRL* RRRNRL RRNNRL 
Q 5/6 RNNNRL* 

7 7/7 RRRRRRN RRRRRNN RRRRNNN RRRNNNN RRNNNNN 
7 7/7 RNNNNNN RRRRNRN RRRNNRN RRNNNRN RRNRNRN 
7 7/7 RNNNNRN RNNRNRN RRNRNNN RNNRRRN 

7 A/7 RRLRLRL NRLRLRL RRLLNRL NRLRLNL NRLLNRL NLNLNRL 
7 4/7 RRLLNLN RRLRLLN N RLLNLN* 

7 5/7 RRLRRLR RRLRRLN RRLNRLN NRLNRLN RRLNLNN* NRLNLNN 
7 6/7 RRRRRRL NRRRRRL NNRRRRL NNNRRRL NNNNRRL NNNNNRL 
7 6/7 RN RRRRL N RNRRRL RRNRRRL NN RNRRL RN RNRRL RRRNRRL 
7 6/7 NNNRNRL NRNRNRL NRRRNRL RRRRNRL NRRNRRL RNNNNRL 
7 6/7 RNNRNRL RRNRNRL 



For the ID circle map, we have determined all allowed periodic sequences up to period 7, 
which are listed in Table fTI|. We have examined their admissibility by using the tangencies 
of the 2D Poincare map and found that ten of these cycles are now forbidden by the tangency 
T5. An asterisk denotes those forbidden sequences in Table fT|. The allowed periodic orbits 
have been located numerically. 
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FIG. 19: Stable period 6 NRLLRN at Q = 1.26055 and the chaotic attractor at Q = 1.26005 in 
the — Xn map. 

VIII. PERIOD WINDOWS 

So far we have discussed only unstable periodic orbits which are embedded in a chaotic 
attractor for fixed parameters. In fact, the symbolic dynamics we have constructed is also 
capable of treating stable periodic orbits appearing while parameters vary once the partition 
lines are determined in the Poincare map. But a more convenient way to determine the 
symbols of a periodic orbit is to use the x^+i — x„ first return map. Take, for instance a 
stable period 6 at Q = 1.26055. We show it and the chaotic attractor at Q = 1.26005 in 



the — Xn map in Fig. |T^. Comparing Fig. |19| with Fig. |14|, one can easily find l=iV, 
4=L, and 5=-R for the two chaotic attractors look very similar. To tell the symbols of the 
rest points we need further to specifically locate the tangencies close to 3 and 6 on partition 
lines •B and •C. As a result, we found a tangency (0.409421208764, 0.513670945215) on 
the left side of point 3 (0.417925742243, 0.515459734623) and the other (0.678580462575, 
0.299669053868) also on the left side of 6 (0.679571155041, 0.299669061327). Hence, the 
letters for points 2, 3, and 6 should be i?, L, and A^, respectively. We have determined 
all stable periodic words with period up to 7 encountered as Q is varied from 0.78058 to 
1.32280 with an increment 0.00001, as listed in Table Naturally, those periods with 
window width less than 0.00001 would be missing. In Table ^ we also list the changes 
of period words with Q, which show the ordering properties: The words in the form ■ ■ ■ C 
would undergo the change from smaller to bigger as Q increases where C corresponds to the 
tangency in the upper branch of attractor on the partition line (see Fig. [1^); while the 
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words in forms ■ ■ ■ and ■ ■ ■ AB would do the opposite, i.e., change from bigger to smaller, 
where is the tangency in the lower branch of attractor along the (see Fig. ^). 



TABLE IV: Period windows with period up to 7 alone the line a = 0.8 + 0.3Q. 



P Range in Q Word and its change with Q 

1 0.78058 - 1.03032 C {R ^ C ^ N) 

2 1.03033 - 1.18422 NC {NN NC ^ NR) 

4 1.18423 - 1.20583 NRNC {NRNR NRNC ^ NRNN) 

6 1.21727 - 1.21783 NRNNNC {NRNNNR NRNNNC NRNNNN) 

7 1.22590 - 1.22596 NRNNNNC {NRNNNNN NRNNNNC NRNNNNR) 

5 1.23171 - 1.23228 NRNNC {NRNNN NRNNC NRNNR) 

7 1.23667 - 1.23673 NRNNRNC {NRNNRNR NRNNRNC NRNNRNN) 

3 1.24492 - 1.24591 NRC {NRN NRC NRR) 

6 1.24592 - 1.24631 NRRNRC {NRRNRR NRRNRC NRRNRN) 

5 1.250853 - 1.25087 NRRNC {NRRNR NRRNC NRRNN) 

7 1.255025 - 1.25503 NRLNNNC {NRLNNNR NRLNNNC NRLNNNN) 
7 1.25645 - 1.256457 NRLRNNC {NRLRNNN NRLRNNC NRLRNNR) 

6 1.26055 - 1.26173 NABLRC {NRLLRN NRLLRC NRLLRR 

NABLRR NLNLRR NLNLRC NLNLRN) 

6 1.26735 - 1.26738 LRRNAB {LRRNLN LRRNAB LRRNRL) 

4 1.27539 - 1.27692 NLNC {NLNR NLNC NLNN) 

7 1.28907 - 1.28909 NLNNNRC {NLNNNRR NLNNNRC NLNNNRN) 
7 1.31289 - 1.31293 NLRNNRC {NLRNNRN NLRNNRC NLRNNRR) 

7 1.31700 - 1.31713 NNRNNLC^ {NNRNNLN NNRNNLC^ NNRNNLR) 
3 1.32043 - 1.32199 LAB {LRL LAB LLN) 

6 1.32200 - 1.32280 LLNLAB {LLNLLN -> LLNLAB LLNLRL) 



IX. PROPULSION ASPECTS 

Apart from its fairly rich non-linear structure, another interesting aspect of the model 
considered here is its use in a mechanical propulsion device. To see these aspects, let us use 
Eqs. (^, ^ to solve for the coordinates (in Cartesian space) of the center of mass of the rod. 
We are particularly interested in this section in the acceleration; two representative graphs 



of the acceleration in the x direction are shown in Figs. ^ and Similar results hold 

for the acceleration in the ?/-direction. These graphs were generated with a particular set 
of parameters and initial conditions of the model whose precise values are not important for 
our purposes here. The interesting aspect is that the acceleration in, say, the x direction 
shows a definite bias in favor of being larger in magnitude in for, in this case, positive values 
over negative values. 

Let us imagine that we place this mechanical model on a cart which is resting on a surface 
and let it begin its motion. The acceleration of the rod will, by Newton's 3'''^ law, cause 
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With friction 

400.0 n 




-500.0 0.0 500.0 1000.0 1500.0 

phi=omega*t 

FIG. 20: Acceleration in the x direction of the center-of-mass of the rod as a function of time (low 
friction) . 



With friction 

270.0 n 




-500.0 0.0 500.0 1000.0 1500.0 

phi=omega*t 



FIG. 21: Acceleration in the x direction of the center-of-mass of the rod as a function of time 
(high friction). 
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a back reaction, causing the cart to oscillate back and forth. First assume the surface is 
frictionless. Even though the acceleration of the rod is larger in one particular direction 
than the opposite one, the cart will not show any net movement - the larger acceleration in, 
say, the +x direction is compensated by the smaller acceleration in the —x directions which 
lasts for comparatively longer periods of time. However, let us now imagine that there is 
some friction between the cart and the surface. It is possible, with the right combination 
of parameters and initial conditions, that this friction "absorbs" the back reaction in the 
direction for which the acceleration is smallest, but doesn't completely absorb the back 
reaction in the other direction. In this case the cart can "push off" the surface, and exhibit 
a net motion. 

Similar effects will hold for motion in the y direction, and so the net motion will be 
a combination of motion in both directions. However, if we place two such mechanical 
devices on the cart and have them rotate in opposite directions, the back reactions in, say, 
the y direction can be arranged to cancel, leaving a net motion in only the x direction. 
Such a device has been constructed (U. S. Patent #4,631,971 and U. S. Patent Application 
7^268,914); one rather striking illustration of its operation is that the device, self-contained 
within a box, when placed in a canoe in a swimming pool, will propel the canoe in one 
particular direction, without the need for a propeller in contact with the water. 

Although perhaps surprising, the principle behind this device is intuitively known to 
many children. Have a child sit on a piece of cardboard on a smooth floor and then tell 
her to move forward by rocking. If one watches closely, the child will rock in one direction 
quickly and return in the other direction slowly. The change in momentum of the child leads 
to a reactive force being transferred to the cardboard, but the uneven rocking results in this 
force overcoming the frictional force in only one particular direction. The net result is that 
the child will propel herself in one direction. 



X. CONCLUSIONS 

In this paper we have studied a simple mechanical model of a periodically driven non- 
linear mechanical system. Bifurcation diagrams, dimensions, power spectra, and Lyapunov 
exponents were derived in order to locate regions of quasiperiodic, periodic and chaotic be- 
havior within the parameter space of the system. Within the parameter space is the the 
coexistence of a chaotic attractor and a limit cycle, which implies that initial conditions may 
determine whether the system will exhibit periodic or chaotic behavior. We also discussed 
some aspects of the model as it relates to a propulsion device. This system, although rela- 
tively simple, exhibits a rich behavior with a number of interesting non-linear effects being 
present. 

We have also performed a symbolic analysis of the model. By constructing the proper 
Poincare section in the phase space for a system of ODEs, the symbolic dynamics can be 
constructed based on the appropriate partitioning of the phase portrait and it turns out 
to be an efficient and powerful way to explore the global properties of the system both in 
the phase and parameter spaces. Up to now, the symbolic dynamics has been applied to 
the analysis of the NMR-laser chaos model [jl2|, [T^, the two-well Duffing equation [jll] 



the forced Brusselator ||^, 0, and the Lorenz model |]T5|, |T6|, [T^, in addition to the system 



discussed in this paper. Along a certain direction in the parameter space this model exhibits 
various properties, such as, periodicity, quasiperiodicity, chaos, ID and 2D features, etc. In 
some other directions or regions of the parameter space the model would also display more or 
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less similar behavior. We have established the 3-letter symbolic dynamics for the model and 
found that the ordering rules of sequences, the forced Brusselator in the regime of annular 
dynamics and the dissipative standard map at some parameters are the same. As a matter of 
fact, the NMR-laser chaos model, the forced Brusselator in the regime of interval dynamics 
and the Henon map with a positive Jacobian also have similar 2-letter symbolic dynamics 
and share the same ordering rules of sequences. It, therefore is meaningful in a sense to 
classify the systems of ODEs according to their ordering rules of sequences. The ODEs 
investigated under the guidance of symbolic dynamics to day are quite limited though. 
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